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摘 要 : 土 体 构 型 对 土壤 水 分 、 溶 质 运 移 过 程 和 作物 成 长 等 有 显著 影响 ,常规 测量 采用 人 工控 土壤 
剖面 取样 、 实 验 室 化 验 分 析 等 方法 ,其 周期 长 效率 低 。 针 对 以 上 问题 ,以 探 地 雷达 波形 及 其 图 像 为 
研究 对 象 , 从 检测 土 体 构 型 的 属性 (层次 、 层 厚 、 土 质 ) 入 手 , 提 出 了 一 种 使 用 探 地 雷达 快速 测量 土 
壤土 体 构 型 的 无 损 探测 方法 。 基 于 探 地 雷达 波形 图 像 的 纵向 梯度 信息 能 够 反应 土壤 分 层 ,采用 包 
络 检 波 法 从 探 地 雷达 回 波 中 提取 包 络 信号 ,利用 Hilbert 分 析 其 瞬时 相位 来 确定 分 层 位 置 ;鉴于 
壤 介 电 和 常数 与 雷达 回 波 振 幅 的 关系 ,采用 探 地 雷达 回 波 振幅 反 演 各 层 介 电 常数 ,由 介 电 常数 推 
雷达 波 在 土壤 中 的 传播 速度 ,以 此 得 到 土壤 剖面 各 层 厚 度 ;根据 探 地 雷达 波形 图 的 图 像 躁 声 与 土 
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壤 砂 壤 比 之 间 存 在 定量 关系 ,提出 采用 主 成 分 分 析 方 法 对 每 一 层 土 质 的 图 像 进 行 噪声 估计 求 得 各 
层 土壤 的 含 砂 量 , 结 合 支持 向 量 机 进而 辨识 各 层 土 壤 士 质 。 建 立 涵盖 地 域 信 息 、 土 壤 指 标 、 探 测 信 


息 、 图 像 多 特 和 


证 。 研 究 表明 


测 思路 。 
关 键 词 : 土壤 ; 土 体 构 型 ; 探 地 雷达 ; 无 损 探 测 


文章 编号 : 1000 - 6060(2022)06 - 1860 - 10(1860 ~ 1869) 


土壤 是 地 球 上 所 有 生态 系统 的 自然 支撑 和 生 
物 屏障 。 近 年 来 随 着 工业 化 ,城市 化 进程 的 加 剧 , 
人 与 土地 之 间 的 矛盾 也 日 着 突出 ,土壤 质量 和 可 持 
续 发 展 已 成 为 社会 经 济 发 展 和 土壤 学 .农学 .环境 
科学 等 研究 的 热点 "。 土 壤 质量 检测 方法 通常 有 安 
观 和 微观 两 种 ,宏观 即 卫星 遥感 探测 法 ,是 利用 卫 
星 遥 感 技术 探测 土壤 表层 属性 ,对 土壤 表层 含水 
K .全 所 和 全 磷 含 量 . 盐 渍 化 .重金属 含 量 等 参量 进 
行 检测 ”” ,能 以 不 同时 空 尺度 不 断 地 观察 地 表情 
况 , 提 供 地表 特 征 信息 * ,适用 于 空间 覆盖 范围 较 广 
的 土壤 表层 信息 监测 ,但 监测 数据 的 空间 分 辨 率 往 
往 偶 低 ;微观 即 点 测 法 ,是 通过 剖面 挖掘 . 打 土 钻 、 


收 稿 日 期 : 2022-03-09; 修订 日 期 : 2022-05-24 
基金 项 目 : 内 蒙古 自治 区 科技 成 果 转 化 项 目 (2019CG071) 资 助 


融合 信息 的 土 体 构 型 知识 库 , 并 编制 快速 识别 土壤 土 体 构 型 的 信息 系统 ,在 内 蒙古 

自治 区 呼和浩特 市 周边 2 个 试验 基地 、6 个 采样 点 .4 类 土壤 土 体 构 型 使 用 该 方法 进行 野外 探测 验 
:在 上 述 地 区 地 表 以 下 1m 范 围 内 的 土 体 构 型 识别 正确 率 达 到 94% 以 上 ;探测 的 每 层 

土壤 厚度 相对 测量 误差 小 于 10%。 该 研究 方法 为 土 体 构 型 快速 检测 应 用 提供 了 一 种 新 的 中 观 探 


取样 方 测 土 样 等 方式 获得 有 限 的 点 测 数据 ”, 来 推 
断 整 个 面 或 三 维 空间 的 土壤 质量 ,测量 准确 ,但 仅 
适用 空间 窗 盖 范围 偏 小 的 检测 。 从 空间 上 来 讲 ,还 
感 探 测 法 与 点 测 法 不 够 全 面 , 咪 需 在 中 观 上 建立 一 
种 客观 ,快速 ,准确 的 土壤 质量 检测 方法 。 

探 地 雷达 是 一 种 利用 高 频 电磁 波 的 反射 探测 
目的 体 即 地 质 现象 的 物探 方法 , 亦 称 地 质 雷 达 , 具 
有 高 效 、 快 速 .连续 无损、 低 成 本 和 高 分 辨 率 成 像 
等 特点 ,已 成 为 公路 路 面 ` 机 场 道 面 .隧道 ,水坝 和 
堤防 等 层 状 结构 基础 工程 设施 无 损 检测 技术 的 重 
要 组 成 部 分 "”。 近 年 来 , 随 着 探 地 雷达 技术 的 快 
速 发 展 逐 渐 应 用 到 土壤 检测 领域 中 ,成 为 宏观 探测 
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与 微观 探测 之 间 尺 度 转化 的 纽带 。 如 吴 志 远 等 
利用 探 地 雷达 早期 信号 振幅 包 络 法 平均 值 方法 和 1 ”研究 区 概况 
时 域 反射 技术 分 别 探测 降水 前 后 野外 农田 表层 的 ee 
BF 人 334 ae AS vy 2 O, H= DA PY Ee N 9 H 
LEBER DARE PE BP IN SIR EHER h I5 EFA ERG 
幅 包 络 值 能 够 获得 与 时 域 反 射 技术 实测 精度 相近 GEEA ei 
IREE SBE A (RDS AC AS OH EE 旗 和 开 2 黄 土 丘 陵 区 的 达 拉 特 旗 ( 图 1)。 土 默 特 左 
F ETIE IRER 等 探 地 雷达 用 于 土 eg So, ee ees 
On eee can 旗 位 于 内 蒙古 自治 区 中 南部 .黄河 北岸 ,属于 典型 
壤 污 染 物 探测 中 ,以 土壤 介 电 常数 为 纽带 ,通过 神 大 陆 性 半 干 早 季风气候。 降水 年 际 变化 大 , 目 年 内 
经 网 络 技术 ,建立 介 电 常 数 、 雷 达 回 波 信号 和 污染 。 分 本 不 区 A 
sE, 7— A 8%, 
土 末 多 个 参数 之 问 的 联系 ,能 够 实现 污染 土壤 中 多 EuR = ae ne 
加 Z^ mmo DY SM AR 
种 污染 物 含量 的 检测 ;也 有 部 分 学 者 将 探 地 雷达 用 en ee eee nee es uaa 
sie OF ede ee Fei (109°54'44.3"E,40°19'46.0"N ) BLUE , Ayie 
于 土壤 深度 和 分 布 的 探测 ,如 赵 艳 玲 等 使 用 探 地 i a 
a cae er cate nse, SSSR, RERI + SRE TE FA RERE 
SHRINE UC BUMBLE ERN LRR 
ZIN Za zB ’ 2 A 
R EARRA EBA AREENA ARER) 砂 nee cena A 
ý px AS = aE +} ha] yey fez aes Ae te 工 JIN] SW o JA TR 
以 上 研究 ,第 一 从 探 地 雷达 回 波 信号 着 手 ,建立 十 。 内 蒙古 自治 区 鄂尔多斯 高 原 北端 ,黄河 中 游 内 蒙古 
二 与 雷达 波形 数据 的 关系 ;第 二 研究 土壤 的 单一 属 a emrt cE EARE 
性 ,如 含水 率 、 层 厚 、 分 布 .污染 物 含量 等 。 本 文 提 一 Lee ues 
a 候 , 受 季风 环流 影响 , 达 拉 特 旗 春 季 降 水 少 ,夏季 降 
出 了 一 种 使 用 探 地 雷达 快速 测量 土壤 
出 条 地 雷达 快速 测量 土壤 土 体 构 型 的 k ue 
da ee ee 水 集中 , 且 容 易 出 现 大 到 暴雨 天 气 ,该 地 区 降水 季节 
无 损 探 测 方法 。 该 方法 分 别 以 探 地 雷达 波形 及 其 。 分配 不 区 AAKE EIIN E 300 
x 2x J; bt AER APE % oh L SY, PAKEM, Y pE eS) mmo 
Se RR TTI aca FRED tit FRERE ET IED 09°53 
法 ,土壤 介 电 常 数 推演 、 基 于 主 成 分 的 图 像 白 噪声 ue 40o161417oN) 附 近 _ 设 3 个 采样 点 采集 的 
PIP RPL EER MEETA aa P NTD ERE E 
与 土壤 物性 之 间 的 关系 ,以 探测 土壤 前 面 的 三 大 主 
要 属性 (层次 JAE .土质 ) ,综合 地 对 土壤 士 体 构 型 2 数据 与 方法 
进行 研判 ,并 编制 了 一 套 快速 识别 土壤 土 体 构 型 的 
信息 系统 。 野 外 探测 试验 表明 该 方法 实现 了 对 试 21 数据 采集 


验 基地 区 域内 土壤 土 体 构 型 快速 地 .无损 地 探测 。 APURE E . 砂 壤 土 / 壤 土 DEWR 
(a) 呼 和 浇 特 市 地 图 (b) 鄂尔多斯 市 地 图 


4 
ir 特 市 y = ip 
Z \ 
: i, 鄂 托 克 旗 Ej 
A ad © F. jg 
d . Ana, BRY 
Lr: i © A j e © a 
i 和 林 格 尔 县 y py & Speen F da 
< ; ias tf > 
SN aes sa F 省 级 行政 中 心 oe F kd we i 
\ 。 Wo 和 和 TS e ce 
{ y P iaaa, 2) ype EE 
01300 km < ce wa AR 0 3120km oe Ean 
a a L_ 一 ~ 河流 湖泊 
a ， 采样 点 


注 :该 图 基于 国家 测绘 地 理 信 息 局 标准 地 图 服务 网 站 下 载 ,其 中 图 a 为 审 图 号 蒙 SC2020)022 号 的 标准 地 图 制作 ,图 b 为 审 图 号 蒙 SC2020)028 号 
的 标准 地 图 制作 , 底 图 边界 均 无 修改 。 

图 1 采样 区 域 图 

Fig. 1 Plot ofthe sampling area 
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体 砂 4 种 具有 明显 差异 的 土 体 构 型 作为 试验 对 象 ， 
在 无 作物 履 盖 的 耕地 开展 雷达 测 线 布 置 和 土壤 剖 
面 点 测 双 盲 测试。 每 种 土 体 构 型 分 别 布置 了 10 条 
长 10 m 的 平行 测 线 。 

点 测 采 样 :在 每 条 测 线 的 中 点 进行 土壤 剖面 控 
掘 ,土壤 剖面 描述 主要 依据 土壤 颜色 .土壤 松紧 度 、 
粒 径 分 布 等 特征 判定 ,测量 各 质地 层次 的 深度 和 厚 
度 , 利 用 质地 层次 排列 顺序 及 层次 厚度 初步 判定 土 
体 构 型 , 取 各 层 土壤 样品 ,进行 实验 室 化 验 测 定 土 
RTEA 

探 地 雷达 采样 :设备 采用 青岛 中 电 众 益 智能 科 
技 发 展 有 限 公司 的 GER-10 型 脉冲 式 探 地 雷达 (图 
2) , 因 试验 区 的 耕作 土 层 厚度 约 为 1 m, 结 合 土壤 介 
电 特 性 受 土壤 质地 影响 的 最 小 范围 ,选取 900 MHz 
的 发 射 天 线 ,雷达 具体 参数 及 指标 见 表 1 。 


图 2 GER-10 型 脉冲 式 探 地 雷达 
Fig. 2 GER-10 pulsed ground penetrating radar 


表 1 GER-10 型 脉冲 式 探 地 雷达 技术 参数 及 指标 
Tab. 1 Technical parameters and indexes of GER-10 
pulsed Ground Penetrating Radar 


技术 参数 指标 


最 大 扫描 速度 215 Traces-s” 

采样 模式 时 间 模 式 、 距 离 模 式 、 点 测 模式 
时 间 窗 范围 1~8000 ns 

增益 范围 -10~160 dB 

中 心 频率 900 MHz 

深度 范围 0.5~1.5 m 

耦合 方式 地 面 耦合 


探 地 雷达 沿 平 行 测 线 往返 探测 数据 ,一 条 测 线 
探测 4 次 。 现 场 测 量 开 始 前 首先 使 用 Pereometer 


V.7 型 介 电 常数 仪 SU-LB 型 土壤 水 份 仪 对 样 地 进行 
探测 ,根据 探测 值 对 雷达 的 采集 参数 进行 设 定 。 当 
遇 到 坑 洼 地 区 时 ,截断 采集 数据 避免 空气 波 对 图 像 的 
影响 。 一 条 测 线 探测 1 次 约 生成 4~6 组 雷达 波形 图 。 
参数 设 定 的 内 容 包 括 时 间 窗 口 大 小 、 扫 描 样 点 数 、 每 
秒 扫 描 数 .A/D 转换 位 数 增益 点 数 等 内 容 ”"。 参 数 设 
置 的 是 否 合理 影响 到 记录 数据 的 质量 ,至 关 重 要 。 

(1) 探测 深度 与 时 窗 长 度 

探测 深度 的 选择 尤为 重要 , 选 得 太 小 会 丢掉 重 
要 数据 , 选 得 太 大 会 降低 垂 向 分 辩 率 。 一 般 选 取 探 
WREE CH ) 为 目标 深度 的 1.5 倍 。 文 中 探测 的 土壤 
深度 为 1 m, 探 测 深度 设置 为 1.5m。 根 据 石和 介 电 
常数 (& ) 确 定 采 样 时 窗 长 度 (Range, 单 位 为 ns), 公 
式 如 下 : 


Hx(e)” 
0.3 

以 达 拉 特 旗 的 砂 土 / 砂 壤 土 为 例 ,该 区 域 为 农耕 
土 ,2 为 11, Range 应 选 40 ns, 时 窗 选 择 略 有 富余 。 

(2) AD 采样 分 辩 率 

雷达 的 AD 转换 有 8 Bit, 16 Bit 24 Bit, H FIR 
测 深度 在 1.5~5m 之 间 ,采集 速 度 平缓 ,选择 16 Bito 

(3) 采集 模式 及 扫描 样 点 数 

采集 模式 分 点 测 模式 .距离 模式 .时 间 模 式 , 试 
验 时 选择 距离 模式 ,该 模式 下 雷达 按照 相等 的 距离 
间隔 采集 雷达 数据 。 扫 描 样 点 数 Samples 有 128、 
256、512、1024、2048 可 选 ,对 于 雷达 中 心 频率 为 
Fa 、 时 和 窗 长 度 为 Range, Samples 应 满足 下 列 关系 : 

Samples > 10-8xRangexFa (2) 

为 保证 高 的 垂 向 分 辨 ,在 容许 的 情况 下 Samples 
尽量 选 的 大 些 。 文 中 雷达 主 频率 Fa 为 900 MHz, 
Range 为 40 ns ,每 一 道 扫描 样 点 数 大 于 360 ,选择 接 
近 值 为 512。 

(4) 扫描 速率 

即 系统 每 秒 采集 的 数据 道 数 ,GER-10 有 32、 
64、128、256、512 共 计 5 种 ,降低 扫描 速率 可 以 提高 
系统 的 信 噪 比 。 当 选择 距离 模式 采集 数据 时 ,扫描 
速率 自动 设置 在 当前 采样 点 数 下 的 最 高 扫描 速 
率 。 如 本 文 主机 连接 900 MHz 天线 ,扫描 样 点 数 设 
置 为 512, 则 扫描 速率 自动 选择 为 512 Traces- s” (É 
每 秒 ), 测 距 轮 的 间隔 设置 为 2 em, 那么 雷达 的 最 大 
工作 速度 为 36.86 km*h', 如 果 天 线 的 移动 速度 超过 
该 速度 时 采集 数据 水 平 距离 的 误差 将 增 大 。 


Range = 2x =6.6xHx (e)” (1) 
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2.2 研究 方法 
2.21 图 像 预 处 理 地 下 所 有 物体 会 产生 回 波 信 
号 ,介质 的 不 均匀 也 会 产生 回 波 信号 局。 土壤 粒 径 
K/N JEKE .含水 程度 .地面 坡 度 等 几何 特性 的 不 
均匀 性 ,使 得 探 地 雷达 接收 的 回 波 信号 包含 多 种 随 
机 噪声 条 波 干扰 。 本 文 利用 均值 法 、 带 通 滤波 
法 .中 值 滤 波 法 .指数 增益 调节 等 信号 处 理 算法 对 
雷达 波 数据 进行 去 除 背 景 影响 ,高 频 杂 波 .来自 地 
面 发 射 的 直达 波 .信号 增强 等 预 处 理 ,提高 信号 的 
信 品 比 "”。 
2.2.2 基于 包 络 检 波 的 土壤 分 层 信 息 识 别 ” 探 地 雷 
达 发 射 的 电磁 波 信号 遇 到 不 同 质地 的 土壤 时 会 产 
生 不 同 波形 幅度 与 相位 的 回 波 ,波形 的 正 负 峰 分 别 
以 白色 和 黑色 标识 .以 灰 阶 或 彩色 标识 。 这 样 , 同 
相 轴 或 等 灰 度 .等 色 线 即 可 形象 的 表征 地 下 反射 
面 。 反射 脉 冲 波形 的 明显 程度 ,是 对 探 地 雷达 图 像 
进行 地 质 解释 的 重要 依据 ”。 

为 抑制 干扰 噪声 的 影响 ,采用 基于 Hilbert 的 包 
络 检 波 法 从 探 地 雷达 回 波 中 提取 包 络 信号 ”", 分 
析 其 瞬时 相位 来 确定 土壤 分 界面 的 位 置 。 探 地 雷 
达 的 回 波 信号 x(t) 进行 Hilbert 变 换 , 公 式 如 下 : 

X(t) =x()* (3) 


RP: XO 为 Hilbert 变 换 后 所 得 结果 ;xO) 为 探 地 雷 
达 回 波 信 号 ;为 时 间 ;* 表 示 卷 积 运算 。 构 建 复数 域 
的 解析 信号 EE) ,公式 如 下 : 

(0) = x(t) +jX(0) (4) 
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AP: e, 为 下 层 介质 的 介 电 常 数 ; esi, 为 上 层 介质 
的 介 电 和 常数; R 为 反射 系数 ;NN 为 表层 土壤 的 总 层 
数 。R 的 计算 公式 如 下 : 

A(t) 


式 中 : 4() 为 反射 波幅 ; A, 为 全 反射 波幅 数值 ,可 以 
通过 金属 板 试 验 获得 。 利 用 探 地 雷达 测量 时 ,顶层 
介质 为 空气 ,第 1 层 介质 为 第 1 层 土壤 。 因 此 可 由 
顶层 的 空气 介 电 常数 ,推演 出 第 1 层 土 壤 介 电 常 数 ， 
依次 可 推出 其 他 土 层 的 介 电 常数 。 而 对 于 介 电 常 
数 为 e, 的 土壤 介质 来 说 ,电磁 波 在 土壤 中 传播 的 速 
度 (w) 可 由 公式 (9) 推 导 : 
pere” (9) 

式 中 :c 为 电磁 波 在 真空 中 的 传播 速度 , 取 3x10"m's 。 
各 士 层 厚 度 计 算 公 式 为 : 

AH. = o; x At) (10) 
式 中 : AH, N38 i 层 土壤 厚度 ; Ae, 为 电磁 波 在 第 i 
层 土壤 中 的 双 程 时 延 。 土 层 厚 度 可 根据 探 地 雷达 
在 土 层 中 回 波 时 延 和 电磁 波 传播 速度 计算 获得 。 
2.2.4 土质 识别 方法 研究 (1) 基于 主 成 分 分 析 的 
土壤 图 像 噪声 估计 。 土 壤 的 粒 径 、 密 度 、 孔 际 度 、 渗 
透 性 .磁性 等 属性 不 同 ,生成 的 回 波 信 号 图 像 包 含 
噪声 不 同 ,因此 通过 估计 探 地 雷达 回流 图 像 噪声 可 
间接 反映 所 测 土质 的 属性 。 陈 晓 东 等 ”对 土壤 的 
探 地 雷达 图 像 进 行 基于 回归 分 析 和 主 成 分 分 析 
(Principal component analysis , PCA ) 的 真实 噪声 均 方 


式 中 :为 虚 部 单位 。 利 用 解析 信和 号 的 实 部 与 虚 部 求 
得 时 域 振幅 包 络 信号 AW) ,公式 如 下 : 


40= [x (+ XO) (5) 
并 求 得 瞬时 相位 信号 o0) ,公式 如 下 : 
d= ein (6) 


x(t) 
然后 对 b(t) 进行 高 阶 最 小 二 乘 拟 合 得 到 拟 合 曲 
线 。 经 反复 试验 , 当 阶 数 取 20 时 拟 合 曲线 的 凹凸 抛 
点 与 土壤 剖面 分 界线 的 实际 位 置 几乎 完全 拟 合 。 
因此 将 该 拟 合 曲线 上 的 凸凹 拐点 处 作为 土壤 分 界 
面 的 时 域 位 置 。 
2.2.3 土壤 介 电 常数 与 土屋 厚度 计算 采用 探 地 雷 
达 回 波 振 幅 反 演 介 电 常 数 方法 “” ,确定 结构 层 的 
介 电 常数 与 回 波 的 波幅 关系 ,公式 如 下 : 
€,=€, [(1+RV0 -RP i=1,2,--,N (7) 


差 的 估计 ,提出 一 种 利用 图 像 噪 声 估计 分 析 土 壤 含 
砂 量 的 方法 ,并 指出 回 波 图 像 噪 声 数值 大 小 与 土壤 
含 砂 量 成 正比 关系 。 本 文 利用 该 方法 ,来 计算 土壤 
的 含 砂 量 。 
探 地 雷达 所 采集 到 的 含 噪 图 像 (/ ) 可 看 作 是 由 
不 含 噪声 的 真实 图 像 (> ) 和 加 性 高 斯 白 噪声 信号 
(n ) 组 成 ,其 表达 式 为 : 
f=rtn (11) 
RP: 为 均值 为 零 TAN 6? WERT AS 
信和 号。 依据 文献 [29 | 的 分 析 , 将 图 像 数据 / 、r 变换 
到 PCA 域 中 得 到 : 
A nC PEA CH (12) 
式 中 :Ai,0 为 求 取 协 方差 矩阵 的 最 小 特征 值 ; C, 
为 f 的 协 方差 矩阵 ; C, Hr 的 协 方差 矩阵 ; e, 为 加 
性 高 斯 白 噪声 信和 号 的 方差 。 由 于 真实 图 像 的 细节 
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信息 不 丰富 且 具 有 信息 元 余 特 性 , r 的 数据 主要 分 
布 在 有 限 主 分 量 上 ,研究 时 可 以 忽略 , 即 存在 如 下 
公式 : 
A,,(C)=0 (13) 
由 此 ,可 以 用 高 斯 噪声 的 方差 来 估算 探 地 雷达 
回 波 图 像 的 方差 值 , 即 公式 (12) 简 化 为 如 下 公式 : 
Ann); (14) 
对 每 一 层 土质 图 像 都 进行 噪声 估计 ,将 计算 的 
方差 值 作为 该 层 土质 识别 的 特征 参数 之 一 ,其 数值 
大 小 可 间接 反映 十 壤 含 砂 量 大 小 。 
(2) 土壤 含水 率 。 在 降水 或 灌溉 与 土 层 深 度 相 
同 的 条 件 下 ,不 同 土质 因 渗 透 性 能 存在 差异 其 土壤 
含水 率 不 同 , 因 此 将 土壤 含水 率 作 为 土质 识别 的 另 
一 特征 参数 。 本 文 使 用 经 典 的 Topp 多 项 式 ™ 来 计 
算 a, 的 第 i 层 土壤 含水 率 ( 9), s; 的 计算 方法 参见 
公式 (7),9 的 计算 公式 如 下 : 
0=-5.3 x 10° +2.92 x 10° x £, -5.5 x 10“ x 
€,+4.3x10°xe? sp) 
(3) 支持 癌 量 机 的 土质 识别 。 以 每 类 土质 图 像 
所 计算 的 6，、9 FSI = Bt RR EI 
CH) CRR EUR PIRRE Ai RAE EE E E 
量 ) 等 组 成 样本 特征 向 量 并 建立 对 应 的 土质 特征 
库 , 将 土质 特征 库 数据 带 入 支持 向 量 机 进行 训练 并 
识别 各 土 层 的 土质 。 


3 BRAD 


3.1 室内 试验 分 析 
在 分 析 研 究 阶 段 ,2019 年 8 月 设计 并 建设 了 土 


(a) 二 维 雷达 图 像 


(b) 数据 处 理 后 图 像 


maraena 


壤 分 析 试 验 仓 。 针 对 土 体 构 型 检测 设计 了 5 层 土壤 
分 层 试验 ,第 1~3 层 和 第 5 层 分 别 铺设 粒 径 为 大 于 2 
cm, F 2 cm、1~2 cm 1l cm 大 小 的 人 砂 壤土 ,第 4 层 
铺设 粒 径 0.5~1 cm 黏 士 层 。 采 集 的 原始 图 像 . 预 处 
理 图 像 与 分 析 结 果 图 像 如 图 3 所 示 , 图 中 纵 轴 为 土 
壤 深 度 , 水 平 纵 轴 为 雷达 扫描 道 数 ,下 文 以 此 图 像 
分 析 土 壤 的 土 体 构 型 。 

将 图 3b 中 每 道 雷达 数据 沿 行径 方向 取 均 值 ,并 
进行 归 一 化 和 去 直流 分 量 处 理 , 得 到 如 图 4 所 示 黑 
色 土 层 轮 廊 线 ,对 其 进行 Hilbert 变换 后 得 到 红色 上 
包 络 线 .绿色 下 包 络 线 。 

图 4 中 包 络 线 的 瞬时 相位 信号 最 小 二 乘 拟 合 后 
得 到 瞬时 相位 拟 合 曲 线 (图 5) ,利用 微 积分 获取 曲 
线 凸 凸 拐点 (图 5 中 绿色 圆圈 ) , 凸 四 拐点 所 在 位 置 
即 为 各 土 层 分 界 点 ,最终 求 得 土壤 分 层 图 像 (图 
3c)。 

表 2 为 室内 土壤 土 体 构 型 试验 数据 结果 。 其 中 
包括 采用 传统 仪器 测量 的 结果 , 简 记 为 测量 值 , 通 
过 本 方法 计算 出 的 结果 简 记 为 计算 值 。 从 表 2 可 以 
看 出 ,本 文 所 述 方法 与 实测 数据 相 比 , 介 电 常数 的 
相对 误差 最 小 为 6.16% ,最 大 为 9.37% , 层 厚 的 相对 
误差 最 小 为 7.58% ,最 大 为 9.13% ;土壤 含水 率 的 相 
对 误差 最 小 为 8.46%, 最 大 为 9.69%。 男 外 , 土 层 的 
雷达 回 波 图 像 噪声 方差 越 大 ,其 对 应 的 土壤 含 砂 量 
越 大 , 且 各 土 层 的 图 像 噪声 方差 差异 较 大 ,土壤 剖 
面 构 型 总 体 识 别 也 较为 准确 。 

3.2 土 体 构 型 知识 库 的 创建 与 系统 搭建 
建立 了 涵盖 土壤 的 地 域 信 息 ( 如 所 属 旗 县 A 
标 范围 .气候 信息 )、 属 性 指标 (土质 、 含 水 率 、 容 重 、 


(c) 土壤 分 层 图 像 


图 3 探 地 雷达 图 像 


Fig.3 Image of ground penetrating radar 
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图 5 瞬时 相位 拟 合 曲 线 及 凸 四 拐点 图 


Fig. 5 Instantaneous phase fitting curve and convex and 


concave inflection point plot 


KERFA R AE) SA fis CR TA 
参数 .探测 时 间 、 探 线 走向 RUE A) 、 图 像 多 特 
征 融合 信息 (噪声 方差 .平均 灰 度 .三 阶 中 心 距 一 


维 炉 等 ) 的 土 体 构 型 特征 知识 库 。 建 立 的 基于 土壤 
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容重 为 1.4~1.6 tm” ,土壤 含水 率 为 4%~7% 的 典型 
沙 地 WIERE REA WEWERE PE 
体 构 型 知识 库 样 ,是 智能 分 析 土 体 构 型 的 基础 
支撑 。 

搭建 基于 探 地 雷达 的 土 体 构 型 快速 检测 系统 ， 
具有 雷达 波形 预 处 理 .数据 分 析 ,智能 检测 和 实验 
报表 自动 生成 等 功能 。 波 形 预 处 理 包 括 取 值 去 直 
达 波 .中 值 滤 波 .高 斯 滤波 .增益 调整 等 ;数据 分 析 
主要 包括 数据 导入 及 筛选 .分 析 检 测 .指标 计算 、 结 
果 呈 现 、 辅 助 功能 ;智能 检测 从 探 地 雷达 数据 中 分 
析 土 壤 的 土 体 构 型 ,并 给 出 各 项 指标 ;实验 报表 包 
括 当前 实验 分 析 数 据 报表 和 历史 分 析 报 表 。 
3.3 野外 试验 分 析 

2020 年 9 月 在 十 默 特 左 旗 和 达 拉 特 旗 采 集 土 壤 
ERA SE EB ED EMIR E EF EF SP OUT 
试验 。 图 6a 与 图 6b 为 壤 士 /黏土 的 分 层 结果 图 像 和 
剖面 图 ,图 6c 5 Al 6d 为 砂 土 / 砂 壤土 的 分 层 结果 网 
像 和 剂 面 图 。 试 验证 明 ,本 方法 计算 的 分 层 层 数 、 
土质 与 实际 测试 的 结果 相同 ,壤土 / 黏 士 的 层 间 厚度 
误差 在 2~6 cm 之 间 , 砂 土 / 砂 壤 土 的 层 间 层 厚 误差 
在 1~6 em 之 间 ,被 检测 土壤 的 土 体 构 型 与 实际 完全 
相符 ,准确 率 约 为 94%。 具 体 的 土 体 构 型 试验 结 
见 表 3。 


4 讨论 


采用 本 方法 野外 探测 的 4 种 土 体 构 型 ,土壤 层 
次 均 为 4 层 , 除 障 碍 层 外 每 层 土壤 厚度 大 约 20~30 
cm, 相 邻 层 的 土壤 质地 具有 相似 性 ,土壤 层次 探测 


表 2 土壤 属性 检测 结果 
Tab. 2 Results of soil property testing 


所 测 参数 第 1 层 土壤 第 2 层 土壤 第 3 层 土壤 第 4 层 土壤 第 5 层 土壤 
介 电 常数 测量 值 /m 3.41 2.86 3.15 3.45 3.44 
计算 值 /m 3.62 3.11 3.43 3.78 3.75 
相对 误差 /% 6.16 8.74 8.89 9.57 9.01 
层 厚 测量 值 /em 20.00 12.00 8.00 45.00 20.00 
计算 值 /em 21.56 12.91 8.73 48.63 21.65 
相对 误差 /%6 7.80 7.58 9.13 8.07 8.25 
土壤 含水 率 Mel% 4.02 2.58 3.42 4.23 4.15 
计算 值 /% 4.36 2.81 3.73 4.64 4.55 
相对 误差 /% 8.46 8.91 9.06 9.69 9.64 
土壤 含 砂 量 。 计算 值 /% 28.27 26.35 25.14 9.25 23.21 
图 像 噪声 方差 值 822.00 792.00 760.00 306.00 680.00 
土质 识别 结果 砂 壤土 WIRE 砂 壤土 黏土 砂 壤土 
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(a) 壤土 /黏土 分 层 结果 图 像 
距离 /m 


时 间 标 尺 /ns 


(0) 砂 土 / 砂 壤土 分 层 结 果 图 像 
距离 /m 


时 间 标 尺 /ns 


深度 标尺 /m 
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注 :图 像 中 的 红色 直线 表示 土壤 的 分 层 位 置 , 红 白 标尺 一 个 颜色 的 长 度 为 10 cm。 
6 探测 结果 对 比 图 像 


Fig.6 Comparison image of probe results 


准确 ;由 于 探 地 雷达 采用 的 介 电 常数 为 仪 融 探测 均 
值 ,输出 的 探测 深度 略 有 偏差 , 相 比 人 工 检 测 土 层 
厚度 误差 在 6 cm 以 内 ;每 层 土壤 图 像 的 噪声 佑 计 差 
异 明显 ,土质 检测 结果 与 人 工 基 本 一 致 ,总 体 土 体 
构 型 识别 较为 准确 。 


本 研究 从 室内 控制 实验 出 发 ,分 析 土 壤 物 性 与 
探 地 雷达 波形 及 其 图 像 间 定性 与 定量 的 关系 ,提出 
利用 探 地 雷达 波形 的 纵向 梯度 信息 分 析 土 层 层次 ， 
以 雷达 回 波 振 幅 反 演 土 壤 介 电 常 数 进而 推 到 每 层 
土壤 层 厚 ,以 图 像 噪声 方差 来 估计 土壤 含 砂 量 的 理 
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表 3 土 体 构 型 试验 结果 表 
Tab.3 Soil configuration test results 


of Tumed Left Banner 


区 域 取得 数据 方式 ” 层 数 ” 层 号 Jam 土质 类 型 
土 默 特 左 旗 ”人工 采样 4 1 30 WE 
2 24 WE 
3 25 WEE 
4 >69 WŁ 
系统 分 析 4 1 36 WŁ 
2 2 «w+ 
3 19 WI 
4 >75 WIE 
达 拉 特 旗 ATR 4 1 9 壤土 
2 22 ”壤土 
3 30 AEE 
4 >70 E 
系统 分 析 4 1 8 壤土 
2 27 ”壤土 
3 24 Ait 
4 >65 FE 


论 ,然后 以 野外 双 盲 探测 的 方法 验证 了 该 方法 在 土 
壤土 体 构 型 检测 中 的 适应 性 。 第 1 次 提出 把 探 地 雷 
达 波 形 图 像 作为 研究 内 容 , 从 图 像 分 析 的 角度 来 识 
别 土壤 质地 ,丰富 并 拓展 了 此 类 研究 项 目的 数据 源 
和 解决 问题 的 思路 ;创新 性 地 提出 将 传统 人 工 检测 
仪表 与 探 地 雷达 相 结合 进行 土壤 物性 探测 ,以 土壤 
的 层次 .厚度 .土质 为 切入 点 综合 地 分 析 土 体 构 型 
的 工程 思想 。 与 以 往 发 表 的 论著 不 同 , 本 文 结合 探 
地 雷达 波形 图 像 分 析 方 法 ,从 综合 的 角度 给 出 了 一 
种 检测 土壤 土 体 构 型 的 方法 ,并 通过 室内 控制 实验 
和 野外 实验 进行 了 验证 ,具有 较 大 的 工程 意义 ; 相 
比 人 工控 剖面 分 析 方 法 , 既 不 损坏 原 有 土壤 结构 ， 
又 极 大 地 提高 了 检测 速度 和 范围 。 

研究 发 现 ,由 于 使 用 的 探 地 雷达 收发 信号 是 脉 
冲 电磁 信号 ,雷达 回 波 数据 受 地 表 植 被 纱 盖 物 .地 
下 植物 根系 等 因素 影响 较 大 ,其 距离 地 面 的 真空 距 
离 大 于 一 定 范围 时 ,将 影响 采集 的 信号 精度 ,为 此 
实际 探测 时 应 避 开 植被 覆盖 时 间或 地 点 ; 探 地 雷达 
回 波 信号 受 水 的 干扰 很 大 ,特别 是 处 于 水 饱和 状态 
时 土壤 的 其 他 特性 在 图 像 中 均 被 掩盖 ,野外 实际 测 
量 应 避 开 雨天 并 在 地 表 比 较 干 燥 时 进行 测量 ,参考 
地 表土 壤 和 角度 ,水 分 等 情况 合理 设置 雷达 探测 参 
数 。 同 时 ,发 现 土壤 障 但 层 含水 率 相 对 较 高 ,依据 
探 地 雷达 回 波 信号 与 水 的 关系 ,可 以 用 这 一 特性 探 
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测 土壤 障碍 层 或 者 水 管 管道 漏水 点 检测 等 。 

本 文 已 经 证 明 在 典型 区 域 使 用 此 方法 可 以 快 
速 无损、 准 确 地 探测 土壤 士 体 构 型 ,但 我 国 区 域 辽 
阔 ,土壤 士 体 构 型 差异 较 大 ,雷达 发 射 信号 在 这 些 
土壤 中 的 穿 透 能 力 、 受 干扰 情况 也 不 同 ,具体 扩展 
到 其 他 区 域 的 试验 需 下 一 步 工 作 中 进行 深入 研究。 


5 结论 


(1) 通过 基于 包 络 检 波 的 土壤 分 层 信息 识别 方 
法 ,可 得 到 不 同 土 体 结构 的 分 层 图 像 , 清 晰 地 展示 
了 地 下 土壤 分 层 特征 ,并 利用 探 地 雷达 回 波 振幅 反 
演 介 电 常 数 、 佑 算 土壤 厚度 , 虽 存在 一 定 测量 误差 ， 
但 土壤 分 层 层 数 、 层 厚 与 现场 剖面 实测 结果 基本 
相符 。 

(2) 实验 结果 表明 :土壤 中 含 砂 量 越 多 ,对 应 探 
地 雷达 波 的 反射 强度 越 大 ,雷达 波 图 像 噪声 越 强 ; 
HRE .培土 . 砂 壤土 . 砂 土 的 含 砂 量 依 次 增多 ,所 
对 应 的 雷达 波 图 像 噪声 逐渐 变 强 。 

(3) 用 经 典 的 Topp 多 项 式 计算 土壤 含水 率 、 基 
F PCA 的 土壤 图 像 噪声 估计 法 来 间接 估算 土壤 舍 
砂 量 多 少 , 最 后 使 用 土壤 含水 率 .图 像 噪声 方差 和 
其 他 图 像 特征 参数 组 成 特征 向 量 , 利 用 支持 向 量 机 
对 土壤 各 层 土质 进行 识别 ,识别 结果 正确 。 
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Nondestructive inspection method for soil profile configuration based on 
ground penetrating radar 


WU Quan', YAO Xijun, CHEN Xiaodong’, ZHAO Min’, ZHAO Huan’, YUN Hao' 
(1. Inner Mongolia Autonomous Region Territorial Space Planning Academy, Hohhot 010000, Inner Mongolia, China; 
2. Inner Mongolia Dazhi Energy Technology Co Ltd., Hohhot 010000, Inner Mongolia, China) 


Abstract: Soil body configuration has a significant impact on soil moisture, solute transport process, and crop 
growth. Generally, measurements are performed using manual digging of the soil profile and sample collection, 
and then taking the samples to the laboratory for testing and analysis, etc., which has a low cycle length and low 
efficiency. To solve the aforementioned problems, by taking the ground penetrating radar (GPR) waveform and 
its image as the research object, and starting from the detection of the level, layer thickness, and soil quality of 
the soil structure, a nondestructive detection method is proposed to rapidly measure the soil body configuration 
using GPR. In this study, the longitudinal gradient information from the GPR waveform image can reflect soil 
stratification by using an envelope detection method to extract the envelope signal from the GPR echo and a 
Hilbert transform to analyze its instantaneous phase and determine the layer position. Given the relationship 
between the soil dielectric constant and the radar echo amplitude, this paper used the echo amplitude of the GPR 
to invert the dielectric constant of each layer, and then to calculate the propagation speed of the radar wave in the 
soil from the dielectric constant, so as to obtain the thickness of each layer of the soil profile. Based on the 
relationship between the image noise of the GPR waveform map and the soil sand to soil ratio, the principal 
component analysis method is used to evaluate the image noise for each layer of soil and obtain the sand content 
of each layer of soil, and combined with the support vector machine to identify the soil quality of each layer. The 
paper establish a knowledge base of soil structure covering regional information, soil indexes, detection 
information, and image multi-feature fusion information, and compile an information system for the rapid 
identification of soil and soil configuration, and use this method for field detection and verification in two test 
bases, six sampling points, and four types of soil and soil configurations around Hohhot City, Inner Mongolia 
Autonomous Region, China. Studies have shown that: (1) the correct rate of soil configuration recognition within 
1 m below the surface of the above areas has reached more than 94%; (2) the relative measurement error of each 
layer of soil thickness detected is less than 10%. This research method provides a new mesoscopic detection 
concept for the rapid detection of soil configuration. 


Key words: soil; soil body configuration; ground penetrating radar; nondestructive inspection 


